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Abstract 

Electronic states near a square Fermi surface are mapped onto quantum 
chains. Using boson-fermion duality on the chains, the bosonic part of the 
interaction is isolated and diagonalized. These interactions destroy Fermi 
liquid behavior. Non-boson interactions are also generated by this mapping, 
and give rise to a new perturbation theory about the boson problem. A case 
with strong repulsions between parallel faces is studied and solved. There 
is spin-charge separation and the square Fermi surface remains square under 
doping. At half-filling, there is a charge gap and insulating behavior together 
with gapless spin excitations. This mapping appears to be a general tool 
for understanding the properties of interacting electrons on a square Fermi 
surface. 

I. INTRODUCTION 

The problem of interacting electrons in two space dimensions is thought to be central to 
an understanding of high temperature superconductivity^ Ordering in three dimensions, 
according to this scenario, arises from the coupling of highly correlated layersi. Attention 
naturally focusses on the properties of the layers and the possibility of unconventional be- 
havior in two dimensions. 

A standard model often studied in this context is the Hubbard model. This begins with 
electrons which may tunnel between nearest neighbor sites on a square lattice, with ampli- 
tude t, and adds an on-site repulsion of strength U. In comparison with the corresponding 
problem in one dimension!, relatively little is known about this model. 



Perturbation theory in U shows very interesting and complicated behavior. With the start- 
ing point of an average occupation number of one electron per site, serious divergences are 
encounteredB. These have been analyzed for more general interactions than the Hubbard 
model interaction^ and are worse than the usual logarithmic type found in other problems. 
Evidently, there is a need for new techniques to be developed. 

The starting point of one electron per site has a square Fermi surface. This paper dis- 
cusses a novel mapping of the square Fermi surface onto one dimensional chains. Two sets 
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of chains arise, one for each axis of the square. Interactions can be decomposed into couplings 
between parallel faces of the square and adjacent faces. These, in turn, result in interactions 
between the chains. Boson-fermion duality§ can be used on each chain separately to find 
the equivalent boson problem, giving a precise definition of this duality for the square Fermi 
surface. 



In a natural way, the interactions separate into purely bosonic terms plus other new op- 
erators. The bosonic terms can be diagonalized exactly and produce behavior similar to 
the Luttinger model in the sense that there are no fermion quasi-particles or Fermi liquid 
behaviorl. There is complete separation of the charge and spin degrees of freedom, as previ- 
ously found in purely one dimensional models!. As a consequence, charge excitations have 
no spin and pure spin excitations no charge. 

Non-bosonic terms produced by this mapping give rise to a new starting point for per- 
turbation theory. This begins with the free electron kinetic energy plus the bosonic terms 
taken together as the unperturbed Hamiltonian. The non-boson terms arising in the map- 
ping are then taken as the perturbation. It is then important to analyze the operators in 
this new perturbation theory, to see if they are relevant or irrelevant, with respect to the 
bosonic starting point. 

This technique of separating interactions into a bosonic piece plus other operators is familiar 
from many other problems, such as the Kondo model!, the backward scattering model!, or 
the sine-Gordon problem in 1D@. It is slightly more complicated here due to the mapping 
of the square Fermi surface onto chains. 

One of the complicating features of the square Fermi surface concerns the modulation of 
the Fermi velocity along the face — in fact it vanishes at the corners. In this new per- 
turbation analysis, it turns out that the operator causing modulation is irrelevant for any 
interaction. The large distance behavior is characterized by a constant velocity. 

It might be argued that doping would destroy these results, since the Fermi surface should 
then depart from a square. That is only correct for nearly free electrons. For sufficiently 
strong interactions, the operator which causes deviations from the square becomes irrelevant. 
A square remains a square under doping, only the value of the Fermi momentum changes. 

Generally, it is found that two different types of operators are relevant for repulsive in- 
teractions. These correspond to Umklapp terms and an antiferromagnetic spin-flip process. 
For a restricted case of repulsive interactions between parallel faces, these extra terms can 
also be diagonalized. They lead to a gap in the charge excitation spectrum and, therefore, 
insulating behavior. There are gapless spin excitations. This range of coupling constants 
also exhibits spin-charge separation. 

There have been studies of coupled ID chains using models for the coupling^ which dif- 
fer from that encountered in this paper. Consequently, the results cannot be applied here. 
Other studies of the square Fermi surface problem have concentrated on models with boson- 
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like interactions!! 2 ! and found non-Fermi liquid behavioilii It is not clear how to extend 
these to include the non-boson interactions which are found to be relevant for short range 
interactions. 

The techniques introduced in this paper are generally applicable for understanding the rel- 
evant degrees of freedom for interacting electrons on a square Fermi surface. Although the 
special case of interactions between parallel faces is emphasized, it should be easy to general- 
ize these results to other models, including the Hubbard model. Future work will doubtlessly 
explore the fall possibilities of this mapping. Some interesting directions are indicated in 
the discussion section. 



II. MAPPING THE SQUARE FERMI SURFACE ONTO CHAINS 

The general problem of interacting electrons in two space dimensions, described by the tight 
binding approximation on square lattices, simplifies if attention is focused on low energy 
excitations near the Fermi surface. For the half-filled hand, the Fermi surface is a square, 
and the important excitations preserve this symmetry. Intuitively, excitations perpendicular 
to a face of the square can be thought of as one- dimensional in character if the momentum 
parallel to the face is zero. For large parallel momentum, this picture is not obvious since 
the Fermi velocity can vary with position along the face. In addition, low energy excitations 
between different faces can also exist, confusing a separation into excitations of each separate 
face. 

These different excitations are contained in the kinetic energy operator of the tight-binding 
model, 

H Q = ~t^2 a\a a jv + ^Yl a\a a icT (2.1) 

(ij)cr ia 

where is the Fermi operator of an electron at site i with spin a , t is the hopping matrix 
element, and the sum over i and j is restricted to nearest neighbors on a square lattice. 
After a Fourier transformation, and rotation by 45°, this can be written as: 

H D = -4t^ cos(k x s') cos(k y s') at^ a^ a (2.2) 

ka 

where s' = s/y2. The fermion field operator is defined before rotation as N^ 1 ap a e tk '' R \ 

for a lattice of N 2 sites, area L 2 , lattice constant "s", and k' = 2TiL~ 1 {n' 2 ,n' y ). The sum is 
the over over the Brillouin zone, and the chemical potential /i, in Eq. fl2.1|) is chosen such 
that the energy of Eq. ( |2.2| ) vanishes for k x or k y equal to ± kp, with kp = 7r(\/2s) -1 . This 
defines the square Fermi surface. 

Consider now that face for k x near kp. The energy can be linearized to find 

Ei(k x ) — v F (k y ) (k x — k F ) (2.3) 
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where vp{k y ) = (2y/2ts) cos(k y s') and the "1" subscript refers to the particular face with 
k x = kp, —kp < k y < kp, and positive (or zero) Fermi velocity. There is a similar lineariza- 
tion for the other three faces, labelled by "2", with k x = —kp, "3" with k y = kp, and "4" 
with k y = —kp, in all cases with positive (or zero) Fermi velocity. 

For a given spin, the Fermi operator describing the "1" face is given, after rotation, by 

ti(x,y) = ± £ aVt e ik vV + ik - x (2.4) 

— * — * 

The k momenta are 45° rotations of the k 1 , k y runs from — kp to kp, x and y are rotated 
coordinates. With R! = s(m,n), these are given by R = (x,y) with \[2x = s(m + n) and 
\ply = s(n — to). The equal time correlation function (tpi[x, y)ip[) is given by 

($(x,vWi) = s 2 I"" e~ ik yV I** ^ e -ik x x-a\k x -k F \ 

J-k F 2lT J-oo 2tt 

(2.5) 

i / ikpx \ 

= irktgU-V) n l f \ ■ \ ) 
F yv ; \2ti(x + ia) J 

where g(l) = 2~ 1 <5/ + (ttI)^ 1 sin(7rZ2 -1 ), with I = n — to representing a discrete "y" coor- 
dinate, and a is a cutoff of order kp 1 for states far away from the Fermi surface. In these 
formula the sums over k have been replaced by integrals, using the invariance of dk x dk y 
under rotation. 



This result is a product of a function of the y variable times a purely ID expectation 
value. It suggests that the Fermi operator describing excitations near the "1" face can be 
written in the form: 

[— 27V 

ti(x,v) = ^ £<?(Z-OVv(z) (2.6) 

where the prime denotes the operator for the square Fermi surface, the operator on the right 
of Eq. ( |2.6| ) describes fermions on independent chains, and \[2y = si. The integer / runs 
from 1 to 2N, since I = n — m with to and n running from 1 to N. 



The expectation value for independent chains 
{w[i 1 {x)'ipxi 2 )) with normalization ipi(x) = 

Ek x ai(kx) e ikxX , is: 



k<kp 

(2.7) 

/ IP LKpX ^ 

$h,l 



2tt(x + ia) 
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which gives Eq. (|2.5|) when £)j g(l — h)g{h — I') = g(l — is used with Eq. ( [2.6| ). 



It is also interesting to consider the mapping of Eq. ( |2.6| ) in the language of Fourier trans- 
forms. For 2N independent chains, the transformed operator is 

27V 

ij> lny {x) = (2N)- 1 / 2 Y / M^) j pyl (2-8) 
1=1 

with p y = 2nn y (2N)~ l and n y in the interval —N to +iV. This corresponds to the in- 
terval in k y from —2kp to +2/cp, since 2kpy = ttI. The discretization of this interval is 
consistent with the spacing of lines of equal k y , which is 2tcL~ 1 2~^, and the phase factor 
k y y = n y (2nL~ 1 2~2)y = p y l. The point with k y = kp has n y = y, and the square Fermi 
surface thus corresponds to 2N independent ID chains with the Fourier transform momen- 
tum restricted to the half-interval. 

An obvious question with this mapping onto independent chains concerns the continuum 
limit used in the normalization of Eq. ( |2.7| ), contrasted to that of Eq. ( |2.4| ). They differ 
by 's/s, which is to be expected for a purely ID system0. For the discrete lattice system 
in ID, the free fermion expectation value (■?/>} (m)-?/>i(o)) behaves as m _1 for large m, while 
{tlj\{x)il)i(o)) for the continuum model behaves as x" 1 . Since x — > ms, the operators must 
differ by the multiplicative factor of ^/s. With the convention used in Eq. ( [2.7] ), the corre- 
lation functions are correctly reproduced by Eq. ( [2.6D , and that is used in this paper. 

This continuum limit cannot be applied to the other dimension, since the Fermi surface 
degeneracy occurs over a finite range in momentum space. The continuum limit can only 
be taken for states perpendicular to the Fermi surface, while the other direction remains 
discrete. 

Similar considerations apply to the construction of the Fermi operator for k x near — kp, 
the "2" states. The corresponding expectation value is found to be: 

t (X; ) = g (l _ (2 .9) 

2tt(— x + la) 

and the ip' 2 field is represented by 

27V 

i;' 2 (x,y)=7ik F 1 Y,9{l-l')^2v{x) (2.10) 
i'=i 

where ip2i{x) is the operator for the "2" fermion on chain /. 

The operator content of this mapping is clear. With the Fourier transform operator on 
the Fermi surface defined by 

V^O) = ipinyix), -f < n y < f 

and: 
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2N 



if> lny (x) = (2N)~^M^ e~ ip » l —N < n y < N (2.11) 
1=1 

It is seen that the Fermi surface operator is just the chain operator restricted to half the 
transverse Brillouin zone. This mapping is useful because interactions between the various 
faces of the Fermi surface have simple representations when expressed in terms of the chain 
picture. This is discussed in the following section. 

Coupling between the chains could arise from the dependence of the Fermi velocity on 
the momentum along the face. This coupling is contained in the k y Fourier transform of Eq. 
and Eq. 



N/2 



Hl= I Y MPyX y 0%V4n» (2-12) 

n y =- — 

Since the sum in this equation is restricted to the half interval, the ipi operators are equal to 



the unprimed operators, and the inverse Fourier transform of Eq. (2.4) can be used to find 

*=E / jV^^W (2-13) 

where 



JY 
2 

ipyQ-V) 



t hV = {2N)~ 1 £ v F (p y )e 

In the following section, it is shown that the off-diagonal elements of are irrelevant 
operators when interactions are included. The diagonal element represents an asymptotic 
Fermi velocity, v a , of magnitude 

,„ = ^£) (2 . 14) 

7T 

and %i> -> v 5ijt. 

Keeping only the diagonal element, the kinetic energy is a sum of independent chains. Each 
chain can be separately bosonized, giving the result: 

H 1 = nv J2 dx pu(x)pu(x) (2.15) 
i=i Jo 

where pu(x) = if)ii(x)if)u(x) , and the product is normal ordered. The Fourier transforms of 
pn{x) are: 



J2 I dx e ip y l+ik * x Pll ( 

i-i J° 

(2.16) 

p u {x) = (2NL)- 1 e- lkxX - ipyl pi(n) 
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With these relations, the kinetic energy for face "1" can be written as: 

77-7; 

H ^ = Wi 5>(*0piH9 (2-17) 

ft 

where n y runs from —N to iV and the sums over n x are to be cutoff as is the ID case. 
The commutator of ip'^x), from Eq. Q2.6| ), with the kinetic energy is 

[^ u (x),H 1 ]=iv -^ u (x) (2.18) 

using either Eq. Q2.12 ) or Eq. (p,15|) for Hi, as a result of boson- fermion for the ID system 
on chain I. 

These steps can be repeated for the other faces of the Fermi surface, including spin, giving 
the total kinetic operator (under the assumption that tiy — ► v 5i^) as: 

n,cr,/3 

where the products are normal ordered, a = ±1 is a spin index, and /3 = 1 to 4 corresponding 
to the four faces. 

The commutation relations are given by 

[pia'{n'),pia{n)) = -[p2a(n'),p 2a (n)] = 2Nn x 5 H _ H ,5 (J ^ 

(2.20) 

[P3a'(n'),p 3 a(n)} = - [p^'{n'), PAa{n)\ = 2Nn y 5fi,-n'Sa,cT' 
With all other commutations vanishing. Fermi operators for the various faces are: 



r— 2N 

y y f g 9{l ~ 



%A x ^y) = \ljr z2 9{i-i)^i\pA x ) 

(2.21) 

nr 2N 

%>A x >y) = \hr 12 g(k-k')ii) k , }l3 >,a(y) 

V k F k'=l 

where a is the spin index, (3 — 1, 2 the face index, I is the discrete coordinate corresponding 
to y, \/2y = si, and x the continuum variable on chain I'. The other two faces are rotated 
versions of the above, with /3' — 3, 4 the face index, k the discrete x coordinate, and y the 
continuum variable on chain k'. 



The necessity to construct fermion operators unsymmetrically in x and y follows from the 
degeneracies over finite regions in momentum space, the flat pieces of the Fermi surface. 
The mapping onto independent chains gives the correct multiple point fermion expectation 
values, since it preserves the correct operator structure. This construction therefore gives the 
correct asymptotic behavior of these correlation functions and the dynamics near the Fermi 
surface. It is particularly useful when interactions between the different faces are included 
as described in the following section. It should be emphasized that no approximation has 
been made beyond linearization of the free electron spectrum at the Fermi surface and the 
assumption that tw is diagonal, as discussed in Section IV. 
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III. DENSITY OPERATORS FOR THE SQUARE FERMI SURFACE 



The interactions between electron are often described by the Hubbard model, 

H = H + UY, n i+ rn- (3.1) 

i 

where H Q is given by Eq. (|2.1| ), U is the on-site repulsion between electrons of opposite 
spins, and nj CT the number operator at site i with spin a. The previous section presented 
a mapping of H a onto independent chains, and it is natural to understand the interaction 
term in this representation. Basic to this construction is the assumption that states near 
the Fermi surface are important, since these are correctly reproduced by that mapping. 



In principle, the procedure to follow is straightforward. An electron operator is just the 
sum of the four operators, ij) = ipi + ip2 + ip3 + i>4, from the four faces. Each can be mapped 
into independent chains, with two orthogonal sets of chains for the two axes of the square. 
Since n i+ = if) + , there are 256 terms, but there are also great simplifications due to 
symmetry and operator dimension. 

The first simplification occurs in the density operators themselves. Consider the contri- 
bution to the spin (+1) density from face "1", using the mapping of Eq. fl2.19|) 



t** 1 - h)g(l - / 2 )<+(x)^ 2+ (x) (3.2) 

hh 

The terms with l x = l 2 gives a purely bosonic contribution, since tp\ ll+ (x)ipii 1+ (x) = pu 1+ (x) 
is of the form discussed in Eq. (|2.16p . Defining this boson piece by p' 1+ (x,y), results in: 

Pi+( x >y) = nk^^g 2 ^ ~ h)pu 1+ (x) (3.3) 
h 

and has the Fourier transform, with conventions from Eq. Q2.16Q , given by nkp 1 p' x , (n), with 

p' l+ {n) = f{n y )p l+ {n) (3.4) 
In this equation, f(n y ) is the Fourier transform of g 2 {l), 

/K) = | *V(0 = \ (i - %l) = i (i - Jal) 

for \n y \ < N, (\p y \ < 7r). Using the commutation relations of the p a ,a{n) operators on the 
right hand side of Eq. (3.4), given in Eq. (|2.17 ), the operator algebra for the boson part of 



the Fermi surface density operator can be found. 



It is interesting to study the equal-time density-density correlation function, for free parti- 
cles, and determine how much is contained in the purely bosonic piece of Eq. Q3.4|) . The 
correlation function is found using Eq. (|2.6j), 
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= ^k-/g\l)l-—-\ (3.5) 
\2tt x + iaj 

while the boson part of the density operator contributes: 

Wx,y)A) ^^k-/Y,9\l-lW) (f — j^V (3.6) 

The sum over V gives (27r 2 / 2 ) -1 + (4)^ 1 5; , with g 2 (/) from Eq. ( |2.5| ). The difference arises 
from the restriction to equal chains contained in p[(x,y). In general, contributions form 
unequal chains must be examined, and if relevant, they should be retained. Nonetheless, a 
purely bosonic piece can always be extracted from a density operator. The unequal chain 
terms are analyzed in the next section. 



The boson contribution to the Hubbard model interaction of Eq. ( |3.1| ) can be extracted 
using this procedure. After the notation used in Eq. fl2.2|) , the contribution from the "1" 
face is given by 



n 



= ^E/ 2 KW(^)pi-(-^) 

n 

Note that the combination Us enters here, and this defines the coupling constant for the 
continuum limit, Us. In the solutions discussed below, only the ratio Us/v Q appears, and 



from Eq. ( 2.12 ), the lattice constant cancels, as in the ID Hubbard model continuum limit. 
An interaction which exists only at one site can be generalized to have different interactions 
strengths between different faces. By symmetry, one expects three different interactions, 
equal faces, opposite faces, and adjacent faces. There is also a complete separation of charge 
and spin degrees of freedom in these interactions, giving six coupling constants in all. 

Charge-spin separation follows from choosing linear combination in the form 

v / 2p a (n) = p a+ (n) + p a -(n) 

(3.8) 

V2cr a (n) = p a+ (n) - p a -(n) 
where a is the face index. These operators satisfy 

[a a '(n'),(T a {n)) = [p a > (ft 1 ) , p a (n)] = 6^5^' 2iV(e a ■ n) 
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where e\ = x, €2 = —x, £3 = +y, and £4 = —y, and all other commutators vanish. An 
interaction term of the form in Eq. (|3.?1 ) can be re-written, using 

ft,af3 

= f<*p(n) [Pa(n)pp(-n) - a a {n)a p (-n)} (3.9) 

n,af3 

provided f a p is even in n and symmetric under interchange of a and j3, which is the case 
here. 



The most general interaction allowed by this construction is then 

H b = J^J2 [p'MA a pp' p {-ti) - a' a (n)B a ^(-n)} (3.10) 

n,a0 

with A aa = V\ and B aa = U\ representing scattering on the same face; Ayi = A 2 i = A 34 = 
A43 = Vi and B12 = B 2 \ = -B34 = -B43 = U2 for opposite faces, while all other elements of 
^4q/3 = V3 and B a p = Us, representing scattering from adjacent faces, are equal. In a similar 
fashion, the boson form of the kinetic energy, Eq. ( |2.19| ), can be separated into spin and 
charge operators, using Eq. (|3lf). The result is 



H = j^j- 22 \p a {n)p a {-n) + a a (n)a a (-n\ (3.11) 

na 

and together with Eq. (|3.10|) , provides a complete description of the bosonic degree of free- 
dom for the square Fermi surface. This provides a starting point from which the relevant 
(or irrelevant) operators are to be determined. 

For the Hubbard model, all of the six coupling constants are equal, with magnitude 2Us. 
This equality is not obviously useful for they may renormalize in different ways when rel- 
evant operators, coming from non-boson terms in the interaction, are included. Although 
the boson Hamiltonian H + Hb from Eq. ( |3.10|) and Eq. ( |3.11|) can be diagonalized, the 



problem is simplified somewhat if coupling between adjacent faces is omitted, V3 = f/ 3 = 0, 
the procedure followed below. 

The parameters describing boson interactions between parallel faces can be more intuitively 
understood by inverting the spin-charge separation. In this way, Hb for parallel face inter- 
actions can also be written in the form: 

H B = f{n y )p a(7 {n) [{V x - + (Vi + Ui)5 a - a >) p aa '(-n) 



l 7T E f{n y ) Pla {n) \{V 2 - U 2 )S a y + (V 2 + U 2 )8 <J ^ rT ,] p 2a ,{-n) ( 3 - 12 ) 



NL . 

na,a 



+ (3 — 4 terms) 
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where a = 1,2. The contribution from the "3" and "4" faces is found by replacing f(n y ) by 
f(n x ), pi by p 3 , and p 2 by p 4 . 

All of these are forward scattering processes. The first term has small momentum transfer 
on the same face, while the second describes small momentum transfer on opposite faces. 
Although this representation makes a physical interpretation more evident, solutions are 
more conveniently found using the separation into spin and charge operators. 

It should also be recognized that the generalized model with six coupling constants is not 
the most general, even for the boson degrees of freedom. In order to have these unequal, 
it is necessary to have a momentum dependent interaction, rather than the on-site interac- 
tion of the Hubbard model. As a result, the momentum dependence along the square faces 
would also appear, which changes the momentum dependence contained in the f(n) factor 
appearing in Eq. (|3.4|). The model which retains only Vi,V 2 ,U2 and U 2 , discussed in the 
following section, not only provides insight into the more general problem, but perhaps can 
serve as a phenomenological description as well. 



IV. RELEVANT OPERATORS AND THE BOSON PROBLEM 

The preceding section discussed the construction of boson operators which are contained 
in the fermion density operators, and showed how to extract products of boson operators 
from the Hubbard interaction. In this section, the boson problem is diagonalized, and the 
remaining terms contained in the interaction, are isolated and discussed. 

Central issues concern the value of the single electron states and the possible breakdown of 
Fermi liquid theory. The conditions for spin-charge separation need to be analyzed, along 
with the nature of the excitation spectrum. Further questions involve the consequences of 
the modulation of the Fermi velocity along the face and the stability of the square Fermi 
surface to doping away from half-filling. These problems can all be studied using the boson 
Hamiltonian, from the proceeding section. 



The starting point for the boson problem is H D + Hb, given by Eq. (|3.10| ) and Eq. Q3.11| ), 
which can be written as: 



H + H B = ^ Yl ( f 1 + V 1 ^) Pa {n)p a {-n) + ( 1 - uj^±\ a a (n)a a (n) 

SSLrr' \ TTV n / \ TTV n I 



n,a 



2 f 2 ( n y)pi( n )p2(-n) - — f 2 {n y )a 1 (n)a 2 {-{n) 



1TV n TTV. 



2V2 „2/ \ 1 \ / ->\ f 2i 



H h{nx)Pz{n)pz{-n) / (n x )a 4 (n)a A (-n) 

HV TTV 

Here the sum over n is defined from —N to N for momentum parallel to the face in question, 
with an exponential cutoff perpendicular to the face, as in the ID problem. The function 
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fain) equals f(n y ) for a = 1,2 and f(n x ) for a = 3,4. Note that the coupling constant 
enters as the ratio V\/v etc., and the lattice constant contained in these objects cancels. 
Terms containing a zero in the denominator are to be omitted. They corresponds to zero 
modes, containing no space-time dependence, and must be treated separately. 

Separation into spin and charge operators is explicit, and each sector can be diagonalized 
separately, using canonical transformations such as: 

p'x{n) = cosh rj (n y ) pi (n) — sinh r\{n y )p 2 (n) 

(4.2) 

p' 2 (n) = cosh f] {n y ) p 2 (n) — sinh rj(n y )pi(n) 

generated by p[(n) = e Sl pi(n)e~ Sl with S\ = (2iV) _1 J2n n x lr l( n y)Pi(n)p 2 (— n). A simi- 
lar transformation diagonalizes the spin degrees of freedom for the 1-2 faces, with 5*2 = 
(2iV -1 ) J2h n^ 1 u { n y) a i(p)a 2 (— n)- The transformation for the 3-4 faces are found by inter- 
changing n x and n y , with 1— >3 and 2— >4 for the boson subscripts. 



H D = —f- [la{n)p a {n)p a (-n) + 5 a {n)a a (n)a a (-n)} (4.3) 



Applying these transformations to Eq. (|4. 1|) leads to the diagonal form 

ttv 
NL 

where a = 1, 2, 3, 4. The charge velocity renormalizations are given by: 

,2^ _ „,2/- _ A , V! f2 . A 2 (V 2 p(n y )^ 2 



1 1. a 



7 (n) = Y 2 (n) = 11+ —an y ) - (4.4) 
V nv J \ nv J 

with 7s(n) = j^n) found from Eq. ( |4.4j ) by replacing n y by n x . The spin velocities, S a (n), 
are given by 7 a (n) with the replacements V\ — > —U\ and V 2 — > —U 2 . In general, the spin 
and charge velocities are different, but preserve the x — y symmetry. 

The parameters used in these diagonalization transformations Si to S4 are given by 



e ~"*"" > = — 7v+vi < 4 - 5 ' 

and rj(n x ) is Eq. (|4.5|) , replacing n y by n x . For the spin degrees of freedom, p(n y ) is given 
by Eq. (|4.5|) with V\ and V 2 replaced by —U\ and —U 2 respectively, while u(n x ) replaces n y 
by n x . 

This completes the diagonalization of the boson problem. It is seen that the spectrum 
is linear in momentum perpendicular to each face, and each will contribute a term linear 
in temperature to the specific heat. Interactions renormalize this linear term, as in the 
Luttinger model. 
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These transformations can be used to calculate the fermion correlation function, and de- 
termine the relevant operators. It is helpful to use the boson-fermion duality on each chain 
separately, representing for example the fermion operator of Eq. ( |2.6| ) by 



iPu*(x) = Zi e -Vu*{x)+ik F x 
^j 2la (x) = Z t e ^2ia{x) - ik F x 
where the "phase" fields are defined by 

VaUr (x) = W'XXWn) e-^-W 

H 



(4.6) 



(4.7) 



with a = 1,2, a = ±1, and separation into charge and spin operators has been used. A 
renormalization constant Zi has been introduced, which contains the contributions from zero 
modes corresponding to the n x = term in Eq. ( [4.7]) . For the purely ID easel, Z = (2ira)^^ , 
but for chains, this can also contain a phase factor exp(i<j) ar7 (l)) corresponding to a uniform 
fermion phase on different chains. 

The corresponding operators for the 3-4 faces are defined by 

Ifetod/) = Z k e-^Uy)+*k F y 

(4.8) 

^ k *(y) = Z k e^Uy)-*k F y 

and these "phase" fields are found from Eq. ( |4.7| ) by the replacements 1—^3, 2—^4, and n x 
by n y . For these 3-4 face operators, y is now the continuum valuable, and k represents the 
discrete x- variable. Thus the plane wave phase factor k x x + p y l in Eq. fl4.7|) is replaced by 
k y y + p x k, for the (p a ka(y) fields, for a = 3,4. 



In the following, it will be shown that fermion operator must be paired on equal chains, 
otherwise the corresponding terms are irrelevant. That simplifies the zero mode factors, 
Z, since products of fermion operators on different chains commute, and potential single 
fermion ordering terms drop out. 

The calculation of a fermion correlation function in the boson ground state separates into 
products of operators in the spin and charge sectors, as seen by the example 

(ip\ l+ (x,t)ip 10+ ) = Z 2 (e^+W e -V?io+> e~ lkpX 
= Z 2 e- %k F % A,(x,t)S,(x,t) 

where 



13 



Aj(a:,*) = /e * e * \ (4.9) 

with time evolution, and expectation values using the p part of the Hamiltonian in Eq. Q4.3|) . 
In identical manner, the B term replaces p\ with o"i, and uses the a part of Eq. 



Diagonalizing the Hamiltonian mixes pi and P2 according to Eq. (|4.2|) , such that A is a 
product involving the "1" and "2" operators 

/ iV -1 2~l n x l coshr](n y )p 1 (n J t)e~ ikxX ~ ip v l -N~ l 2~^^n~ 1 cosh r](n y )p 1 (n) 
A h (x,t) = (e n e n 

(4.10) 

with A 2 (x,t) given by replacing coshr](n y ) by sinh.r](n y ) and p\ by p 2 - 

In this expectation value, the time evolution pi(n,t) = p\{n) = e^^M 1 i s to be used. 
These can be calculated using the Baker-Hausdorf formula to find 

A u {x,t) =exp|-(4iV)- 1 J2 E n ~x l cosh2 V(n y ) e~ ak * (1 - e W + ^ - %7(%)*)) I 

[ n y =-N n x =l ) 

(4.11) 

with ^2i(x, t) replacing x by — x and cosh 2 ?]^) by sinh 2 i](n y ) . Converting the sum over 
n x to an integral over k x , it is seen that this integral diverges to (— oo) at k x = 0, unless 
I = 0. Thus, at t — 0: 

1 + dp 

ia \ 2 



A u (a; > t) = «J l , — — (4.12) 
Va; + za/ 

where c? p = iV -1 / dn y cosh 2 ^(n y ) . The time dependence is complicated by the dependence 
of velocity on n y but as long as it remains positive, Eq. (|4.11|) provides a calculation of the 
operator dimension of A\, since t scales as x. Proceeding to the other expectation values, 
there results: 



A 2 i(x) = 6i >0 y- 



d p 



■i a 



x + %a 

B u {x) = Si t o 



l + d a 

ia 



,x + ia 

B 2 i(x) = 5, 



d n 

ia 



—x + ia 



and consequently 
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d p + d a 

%ol \ I a 2 '' 



<*U»)*u»> = ft, {—) {—) (4.13) 

where d a is found as in Eq. (|4.12|) for the spin degrees of freedom replacing rj by \i. This 
is the result of independent chains of the Luttinger model with a suitably chosen coupling 
constant. The time dependence is more complicated since the boson spectrum depends on 
n y , but the essential physics remains ID in character. 

The occupation number for an electron (nn a (k x )) on chain / is an integral^ over the cor- 



relation function of Eq. (|4.13 ). It is seen that only fermion operators on the same chain 
contribute, so the occupation number is independent of I. Since the x dependence of Eq. 
( [4. 13] ) is the same as in the Luttinger model, it can be concluded that (nu a (k x )) is continu- 



ous at k x = hp. For other faces, the result is identical. Bosonic portions of the interactions 
on a square Fermi surface destroy Fermi-liquid behaviors. 

It is now appropriate to check the assumption of Eq. (|2.14f) concerning the off-diagonal 
matrix elements of t^, using the boson problem of Eq. (|4.3|) as the starting point. This is 
done using perturbation theory in t^y to check whether this operator is relevant. One of the 
terms to be studied is 

T la (x) = J2 t lv 4 la {x)^ 1V(J {x) (4.14) 

and in second order perturbation theory would contribute 

Jdx Jdt J dx' J dt'(T la (x, t)T la (x', £')> (4-15) 
where time evolution and averaging is with respect to H + Hb- 



Using the boson representation for the ip a ia{x) operator, the expectation value in Eq. (|4. 14|) 
is calculated to be proportional to: 

, / 2 \d[h-h) 



(R-V'r + ia) 4 \R 2 -V'\ 2 + a 2 



where V is an average (finite) velocity, instead of doing the integral over / -f(n y ),R = x — 
x',r — t — t', and d(l\ — h) is given by: 

f N 

d(h — h) = N~ x j dn y sinh 2 r)(n y ) + sinh 2 /J,(n y ) [1 — cosp y (li — l 2 )] (4.16) 

Jo 

For large l\ — h, the cosine term gives a vanishing contribution, and Eq. ( |4.15| ) reduces to 
the expected result from naive dimensional analysis, using Eq. ( f4.12|) , with d(oo) = d p + d a . 

Since d(l± — l 2 ) > 0, the integrals over space-time in Eq. ( |4.15| ) are convergent at large 
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distance giving a factor of the two dimensional volume. The sum over l\ and I2 is propor- 
tional to N, since U^i 2 falls off as (Zi — k)^ 2 - This result means Tu ia (x) is an irrelevant 
perturbation, the usual result in two space-time dimensions, when the operator dimension 
is greater than two. The reduction from two to one space dimensions is a consequence of 
the finite sum over chains, a result reminiscent of independent one dimensional chains. 

In similar fashion, the contributions from the other operators can be seen to be irrelevant. 
Consequently, the modulation of the Fermi velocity along the face of the Fermi surface is an 
irrelevant perturbation about the boson Hamiltonian@. This result should be contrasted with 
perturbation theory about the free fermion Hamiltonian, which contains (log) 2 singularities 
resulting from the vanishing Fermi velocity at the corners of the square Fermi surface. 

Any non-zero interaction will lead to the irrelevancy of the Fermi velocity modulation, but 
a much more restrictive condition results for the stability of the square Fermi surface itself, 
away from half-filling. The appropriate perturbation is the number operator which is given 
for the "1" face by: 

AfiN la = A/^tt^ 1 ) [dxJ2 9(1 ~ h)g(l - h)4 h a(^iha(x) (4.17) 

l,hh 

where Afj, is the change in chemical potential from half-filling. The corresponding correlation 
function to second order is proportional to: 

1 / a 2 \ d (h~k) 

Yl 9 2 (l - h)g 2 (l - k)-rs — 77, — — -T2 m — t//2 2 ; — 2 
t j^ l2 (R — V't + 1a) 1 \R 2 — V t 1 + a 2 J 

since the two derivative operators are absent. For this to be irrelevant, d(l\ — h) > 2 for 
all l\ — I2, as discussed above. This requires a sufficiently strong interaction to be present. 
Under this circumstance, the only contribution is the diagonal term, h = h, given by 

F h 

where the factor of 1/2 comes from J2i 9 2 {l ~ h)- The presence of this term simply shifts 
the Fermi momentum uniformly, the square Fermi surface remains a square. 

Strong interactions tend to localize single electrons to a chain, at least in lowest order. In 
higher orders, or when interactions between adjacent faces are included, it is possible that 
the situation might change. It seems intuitively plausible that a crossover from a square to 
a rounded Fermi surface should occur, but it is not yet clear which terms are responsible for 
this. Intra-facial hopping terms, such as etc, would lead to rounding. These do not 

appear in lowest order, but could enter in connection with some of the relevant operators 
discussed in the following section. This problem requires further attention. It is nonetheless 
very interesting that the square Fermi surface is stable at this level. 



16 



V. RELEVANT INTERACTIONS 



Interactions between chains are generated through the mapping of the fermion operators 
for the four faces onto independent chains. As described in the previous section, it is pos- 
sible to extract a purely bosonic interaction, and use this to define a new starting point 
for perturbation theory. Other interactions remain, and it is necessary to construct a new 
interaction Hamiltonian, and attempt to solve this new problem. This situation is similar 
to the Kondo problem or the backward scattering model, where a similar reconstruction of 
the interactions is possible, leading to a new expansion. 



Consider the interactions involving parallel faces, the "1" and "2" operators. The lattice 



term t/^ni+n,;. would contain 



J2n i+ n^^J2 I dx */+( a; )*i,+ ( a; )*l-( a; )*i,-( a; ) ( 5 - X ) 
i i J 

with the ^i ;a (x) = ip[i a (x) + ip^ix), and these primed operators are, in turn, sums over 
independent chains given by Eq. (|2.6| ). There are similar terms for the "3" and "4" faces, 
and cross terms containing interactions between adjacent faces. Cross terms are not consid- 
ered here. The task now is to extract the purely boson interaction terms in the sum, giving 
Hb of the previous sections, and find the remaining operators. These new operators then 
provide the desired perturbation expansion. 

A conclusion of the previous section is that the sum over independent chain Fermi oper- 
ators must have chain indices paired. If a single unpaired index appears, that operator will 
contribute an operator dimension greater than 1/2, as in Eq. ( |4.13| ), and the remaining three 



greater than 3/2, giving an overall dimension greater than 2. In the second order correla- 



tion function using Eq. (|5.1|) , the overall dimension will be greater than 4, and hence is an 
irrelevant term. This is another manifestation of the irrelevance of single particle hopping 
for the square Fermi surface. 

Consider first the terms which occur when fermion operators with the same spin are paired 
on the same chain. If they have the same face index, there are boson operators already 
included in Hb- Otherwise, there are two types of terms 

H BS = Un 2 kp 2 j dx J2G(h-l 2 ) W lh Jx)iP 2h + (x)VL_(x)V„ 2 _(z) + h.c.l (5.2) 



hh 



Hu = U^k~ F 2 \ dxJ2G(h - k) [< + (x)^ 1+ (x)4 2 _(x)^ 2 _(x) + h.c] (5.3) 
hh 



where 



G(h - k) = £ <? 2 (Z - h)g 2 (l - k) = {UY l 8h,h + 2- 1 7T- 2 (Z 1 - k) 
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Physically, Hbs describes backward scattering on same or different chains, while Hjj is an 
Umklapp process. 

In a ID system, it is known that Hjj is relevant and Hbs is irrelevant for repulsive 
interactions^!, while these roles are reversed for attractive interactions. The situation is 
similar here, as can be shown using the fermion-boson duality principle from the previous 
section. 



The Umklapp term for li = l 2 in the boson representation is given by: 

Hv = -UG(0)(2k F a)- 2 Y: J dx \ e ~ Alk ^ (M*) + W(*))+hc 



(5.4) 



where the <pi p i operators are defined in Eq. ( |4.6| ). Hbs is the same with the <pi p i operators 
replaced by tpxoi operators, and the e ikpX phase factor is absent. The minus sign results from 
ordering Eq. Q5.3p into the form (— )4 ! \i + 4 ! 2i-' l l ) ii-' l l J 2i+- ^ can then be confirmed that the ma- 
trix element of this operator between the states (ipli (x 2 )ipu 1 +(xz) I an d V , 2/i+( x i)V , i«i-( x i)) 
is the same as the boson equivalent in Eq. (|5.4j ), sandwiched between the boson equivalents 
of these two states. 

Diagonalizing H a + Hb as in section IV with Si, will transform the exponent in Eq. 
according to 

e 51 {2N)- 1 n' 1 (pi(n) + p 2 (n)) e -***-^ e -Si 

n 

= {2NY X Y. n- 1 e-^ ) { P i{n)+p 2 {n)) e ~^-W 

n 

from which the scaling dimension of the Umklapp operator, du is seen to be 

1 N 

du = ^ E (2e- 2 ^ ) ) (5.6) 

n y —N 

as discussed in the previous section. This means the correlation function appearing in sec- 
ond order using Eq. ( |5.3|) falls off as (distance) ~ 2du . The zero interaction case, rj = 0, has 
2du = 4 and for repulsive interactions du < 2 from Eq. (|4.5|) . This term is then relevant 
and must be retained. 



(5.5) 



The phase factor e 4lkFk for half-filling is evaluated to find Ak F x = A(ix2~ 1 l 2 s^ 1 )(s 2 _1 / 2 )(m + 
n) — (mod 2ir) and makes no contribution. Notice that for strong interactions, such 
that the Fermi surface remains square as discussed in the previous section kp shifts away 
from 7rs _1 2~2 under doping and the Umklapp term describes a (quantum) commensurate- 
incommensurate transition!^. This transition will differ somewhat from the one dimensional 
problem due to the parallel momentum dependence contained in Hb- 
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Using the diagonalizing transformation for the spin degrees of freedom contained in Eq. 
(5.2), the operator dimension for equal chain backward scattering is given by 



N 

l BS 



1 N 



2N 

and Eq. (|4.5|) shows dss > 2 for repulsive interactions. For attractive interactions, the roles 
are reversed, djj > 2 and dss < 2. 



When li l 2 in Eq. (|5.3|) , the scaling dimension is changed, becoming dependent on (7i — l 2 ) 
through the k y dependence in Hb- In addition, the spin degrees of freedom do not cancel in 
the exponent when the fermion operators are bosonized and contribute to the exponent of 
the power law. The result for the operator scaling dimension is found to be: 



1 N 

du{h-h) = — £ [ e -^K)(i + e ^-< 2 )) 



N 

2N 

Tly— — JV 

_|_ e -2/i(n H )Q _ e ip v (h-h) (5-8) 



For repulsive interactions, rj > and /i < 0, and the value of d(l\ — l 2 ) is not necessarily 
small. The result for dss{h — h) is found from Eq. ( |5.7| ) by interchanging r\ and \i. 

At this stage of the analysis, it is helpful to reflect on the significance of the transfor- 
mations which change the operator dimensions, as determined by rj and /i. For general 
repulsive interactions, rj > and fi < 0, from Eq. ( |4.5| ). For these signs, and sufficiently 
large magnitudes, dssih — h) > 2 for all h — h, and du(l\ — 1%) > 2 for li ^ l 2 . These terms 
are irrelevant. The surviving relevant operator from this pairing is given by chain diagonal 
Umklapp processes, with dimension dy < 2 from Eq. (|5.6|). 



With the coupling constants chosen to lie in this range, it is possible to anticipate the 
existence of a solution. This will be constructed by dropping the irrelevant operators, and 
expanding the relevant operators in the small parameters, e~ 2v or e 2M . There is a similarity 
here to the ID massive Thirring model, where expanding, the mass term when the analogous 
parameter is small, leads to the Klein-Gordon equationt3. 



It is necessary to examine other pairings of Eq. (|5.1[) to extract further relevant opera- 
tors. Pairings which involve the first and third terms (with second and fourth also paired) 
contain terms like 

V'i t « 1+ ^ 1 -V'2i 2 -V'iz 2 +, iAh+^lh-^ih-Ah+i and ^\ h+ ^\ h -^ 2 i 2 -^ 2 i 2+ - After bosonizing these, 
it is found that they are all irrelevant. The first two give the operator dimension of Eq. Q5.8D , 
but with rj — > — rj and \x — > —fi, and these are larger than two in the coupling constant region 
of interest. The third resembles single particle hopping with an extra v2 in the exponent 
after bosonization, and is irrelevant. 

Cross terms which leave a factor of e 2tkFX after bosonization can be neglected since they 
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have a momentum transfer of 2kp. This would require a large energy transfer and the in- 
termediate states would lie far away from the Fermi surface. 

The remaining pairing, with the first and fourth terms on the same chain together with 
the second and third, contains relevant operators. These are given by — ip\i 1+ ip 2 i 1 _ipli 2 _ip 2 i 2+ 
and -V>n 1+ / 2 Ji-Vi 2 -V'H 2+ > plusthe (1 <-> 2) interchange. Terms such as ipl h+ ip lh _ipl h _ip 2h+ 
are irrelevant, and those involving e 2lkpX after bosonization can be discarded as previously 
discussed. 



Both relevant operators involve products of 2k p spin-flip excitations on different chains. 
(The l\ = l 2 terms has already been counted in the first pairing, characterized by djj , dss 
and a boson term). After bosonizing, these terms give respectively: 

H DU = -U(2k F a)~ 2 f dxYl G(h - h) le^M + <P P h(x) + <P*h(?) - <p ah (x) + (1 <_> 2 ) 

(5.9) 



H DB = -U(2k F a 



-2 



dx J2 G(h-h) \eW Xi 



<P ph {x) + <p all {x) ~ 0ah( x ) + (1^2) 

(5.10) 



Here, \j2ip pl {x) = (fi p i(x) + cp 2 . p i and \p2(p ai = ip la i(x) - <p 2a i(x) with fields <pi a i(x), etc, 
defined in Eq. ( |4.6| ). It should be noted that (p p i(x) involves the combination [pi(n) + p2{n)] 
which is multiplied by e~ r,(jly ^ under diagonalization, as in Eq. ( |5.5|) . Similarly, the (p a i{x) 
involves (o"i(n) — c 2 (n)) which is multiplied by e^ 71 ^ under diagonalization, from Eq. G5.5|) . 
Both multiplicative factors should be small for the corresponding operator to be relevant. 
Note that the fermion (1 <->• 2) interchange requires p p i to change sign, since there is a sign 
difference in the boson representation of ipi an d ip2, while (p a \ does not change sign. Under 
hermitian conjugation both change sign. 



It is appropriate to label the first of these terms as "different chain Umklapp", Hp,jj } while 
the second is a "different chain backward scattering", Hp>B- Their scaling dimensions, after 
diagonalizing Hb, are given by 



12, 



1 iv 

— [e~ 2v(ny) (l + e^ 12 ) 

N n y =~N 

_|_ e 2 M("!/) _ e ¥yh2^j 



(5.11) 



d DB {l 



12, 



1 JV 

^ E [«-*-> (1 

ny=—N 

_|_ e 2 M(™I/) A _ 



JPyll 



(5.12) 
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Both are relevant (dimension less than two) for 77 > and fi < 0, the repulsive case of in- 
terest. There is a region of parameter space where all three relevant operators have a small 
dimension. In this case, the exponents contained in Eq. (|5.4|) , Eq. (|5.9|), and Eq. ( |5.10| ) can 
be expanded. 

Scattering between parallel faces results in an equivalence to interacting fermions on sep- 
arate chains. In this section, the interaction terms have been extracted and their scaling 
dimensions calculated. There is a strong coupling region, characterized by three operators 
with small dimension, d Uy d B u and d^B, which must be included. The following section 
presents the construction of the solution. 



VI. GAP GENERATION AND THE INSULATING PHASE 

This section begins with the boson problem of section III, adding the terms with relevant op- 
erators from the previous section. It would be interesting to find a solution for this problem, 
however this has not been found for general coupling constants. But there is a limit where 
a solution is possible for particular choices of coupling constants. This choice corresponds 
to the parameters e~ 2v ^ ny ^ and e 2 ^™^ small, such that the exponents of Eq. (|5.4j ), Eq. ( |5.9|) 
and Eq. (|5.10| ) can be expanded. Physically, this is a strong repulsion regime, with some 
restrictions on the spin dependent interactions, U\ and U 2 of section III. 



The boson part of the problem, H Q + Hb is given by Eq. (|4.1| ). From the preceding section, 
the three relevant operators must be added to define a new "interaction" term. Diagonaliz- 
ing H a + Hb, using the S operators of Eq. ( |4.2|) , transforms the "phase" fields appearing in 
the relevant operators, generally decreasing them (for e~ v and e M both small). Expanding 
these phase factors to quadratic order gives the new interaction Hamiltonian. For the 1-2 
faces, this is: 



H = H D + H int 

H int = -2U(2k F a)- 2 J dx Y^Gih-^e'^l^ + ^l^ + i^-^^e- 3 t 6 - 1 ) 

where H D is Eq. ( |4.3| ) and S = Si + S 2 from Eq. ( |4.2| ). Not only do the first order terms 
cancel, but the cross terms between spin and charge operators cancel as well. To this order, 
the Hamiltonian separates into independent spin and charge sectors. 

After Fourier transforming, the new Hamiltonian can be written in terms of the density 
operators as H = H p + H a , where: 

H p = 2jy^S\TK) [Pi(n)pi(-n) + p 2 (n)p 2 (-n)] 

ft 

+ M 2 p (n y )k- 2 \p x (n) + p 2 (-n)\ [ Pl (-n) + p 2 (-n)]} (6.2) 
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Ha = ^jrJ2\ 5 ( n y) \vi{n)(Ji{-n) + a 2 (n)a 2 (-n)] 

n 

+ Ml{n y )k~ 2 [a^n) + a 2 (n)) [a^-ft) + a 2 (-n)}} 
with velocity renormalizations given by Eq. (t4.4j), and 



(6.3) 



M 2 (n y ) = e-Mn y ) m 



M 2 a {n y ) 



A-Kv U5(n y ) e 2^{n y ) 
(fci?a) 2 



(6.4) 



/ 2 (o)-/ 2 K) 



Here rj and [i are given by Eq. (|4.5| ) and f 2 (n y ) under Eq. ( |3.6|) . The combinations in these 
M 2 are chosen such that the spectrum of Eq. ( |6.2j ) and Eq. (|6.3| ) takes the "relativistic" 
form 



E 2 (n) = k 2 x v 2 ol 2 (n y ) + M 2 (n y ) 



El{n) = k 2 x v 2 5 2 (n y )+M 2 (n y ) 



(6.5) 



as can easily be confirmed by diagonalizing these two. The contributions from the 3-4 faces 
result in an identical Hamiltonian, with r] x «-> n y . 



It is clear that the charge density gap is only weakly dependent on the transverse mo- 
mentum, through the velocity renormalization. On the other hand, the spin density gap 
squared vanishes aX n y = 0, and increases linearly with momentum. A quadratic depen- 
dence on n y would have been intuitively understandable, since it could be renormalized to 
give an isotropic excitation spectrum, proportional to n 2 + n 2 . Perhaps the inclusion of 
interactions between adjacent faces could change \n y \ to n 2 , but that remains unclear at 
present. 

Renormalization of the bare coupling constant U in Eq. (|6.1| ) is required, if the limit a — » 
is to be taken, as in the purely ID problem!. Alternatively, one can regard the M 2 quan- 
tities, as physically observable, setting a length scale for which the bare parameters are 
determined. 



This completes the present analysis of interacting electrons on a square Fermi surface. It 
focussed on interactions between parallel faces, although the method can be readily general- 
ized to adjacent face interactions. Perhaps the most interesting result concerns the mapping 
onto chains, and the ability to extract relevant operators in the chain picture. 

It should be possible to study this mapping more systematically using renormalization group 
scaling. If each non-boson operator is assigned a coupling constant, as in ID systemsi, these 
will satisfy scaling equations. Presumably, some will scale to stronger coupling strengths, 
and solutions such as those presented here, could characterize the asymptotic behavior. It 
will be interesting to confirm, or disprove, the strong coupling region outlined here. 
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Doping away from this point will generate a phase factor in the Umklapp terms, Hy and 
H DU , but not in H DB . This phase factor of exp^iSkpx), where Skp is the departure from 
the half-filled value, changes the nature of the problem. The low energy excitations then in- 
clude some charge states, and it might be expected that a new type of generalized Luttinger 
model results. 

Spin excitations are more problematic in this solution. Although the spin excitation energy 
vanishes when p y — > 0, it does not vanish linearly. At k x = 0, the spectrum is proportional 
to IPyl 1 / 2 . It is not clear if this is an artifact of parallel face interactions and will be modified 
by adjacent face scattering. Another possibility is contained in the problem of uniform phase 
ordering in the chain mapping. As indicated under Eq. ( [4.7| ), a uniform phase on each chain 
is permitted. Perhaps other configurations can be found with lower ground state energy 
than the present result. 

VII. DISCUSSION AND CONCLUSION 

This paper has presented a mapping of the square Fermi surface onto two sets of chains, one 
set for each axis of the square. It correctly describes electronic excitations near the half-filled 
band, with the average occupation number near one. Interactions between electrons sepa- 
rate naturally into boson operators on these chains, plus additional operators. The kinetic 
energy plus the boson operators define a new starting point for perturbation theory, with 
the additional operators as the perturbation. 

This method provides a very powerful new tool for analyzing the properties of the square 
Fermi surface. The purely bosonic part of the problem can be diagonalized exactly, giving 
behavior similar to the Luttinger model. It exhibits power low correlations functions and 
has no jump in the occupation number at the Fermi level. Modulation of the Fermi velocity 
along the face of the square is irrelevant. Rounding of the square is similarly irrelevant for 
sufficiently strong interaction strengths. 

All of these results are derived using the new perturbation theory about the boson part 
of the problem. The interactions also contain non-boson terms, which have been analyzed 
using boson-fermion duality for the chain representations. There is a region of strong cou- 
pling when the relevant operator generated by interactions can be solved exactly, in the 
sense that the sine- Gordon problem becomes the Klein- Gordon problem for strong fermion 
interactions^. For half-filling, a gap appears in the charge spectrum, while the spin excita- 
tions remain gapless. 

While the mapping works for very general interactions, this paper has concentrated on 
a special case of parallel face interactions. For a complete solution, interactions between 
adjacent faces and possible momentum dependences need to be included. A convincing so- 
lution of the Hubbard model awaits an analysis of these generalizations. 
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When doped away from half-filling, spinless fermion states above the charge gap are oc- 
cupied. This situation is reminiscent of the commensnrate-incommensnrate tranaitioni, 
with additional complications due to the transverse momentum dependence. Although the 
spin degrees of freedom are unchanged, it will be interesting to explore the changes in the 
charge excitations. In this way, the very interesting issue of superconducting correlations 
can be resolved. 

Lattice degrees of freedom must also be included, and a logical procedure would involve 
diagonalizing the interesting electron problem first. Adding in the electron-phonon cou- 
pling, selecting out the relevant operators, and solving this new problem, would cast new 
light on the problem of superconductivity in the Cuprate compounds. It seems these steps 
are all possible within the mapping framework. 

It would also seem this problem of interacting fermions on a square Fermi surface could 
be addressed in other ways. The requirement that fermions must be paired on separate 
chains could be realized by a gauge symmetry. Since the single chain problem can be solved 
by a Bethe AnsatzEl, perhaps this new model of coupled chains has a solution involving 
superpositions of these single chain Bethe Ansatz hypothesis wavefunctions. It might be 
simpler to test this on the spinless fermion problem. 

Solutions for the problem of interacting electrons on a square Fermi surface, of the type 
presented here, represent a new regime with correlation effects dominant. Much of the in- 
tuition based on Fermi liquid behavior, is no longer applicable. If it can be shown that 
the high temperature superconductors fall in this regime, this is truly strongly correlated 
superconductivity. 
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